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Abstract In each of the 10 cases with propagators of unit or zero mass, the finite part 
of the scalar 3-loop tetrahedral vacuum diagram is reduced to 4-letter words in the 7- 
letter alphabet of the 1-forms Q := dz/z and u p : = dz/(\~ p — z), where A is the sixth 
root of unity. Three diagrams yield only ((Q 3 uo) = ^ti" 4 - In two cases 7r 4 combines with 
the Euler-Zagier sum ((fl 2 u 3 u ) = J2 m >n>o(~ ^) m+n /m 3 n; in three cases it combines with 
the square of Clausen's Cl 2 (vr/3) = Q : C(^ a; i) = Y,n>o sin(7m/3)/n 2 . The case with 6 
masses involves no further constant; with 5 masses a Deligne-Euler-Zagier sum appears: 
3? ((Q 2 U3UJi) = J2 m >n>o(~ l) m cos(27m/3)/m 3 n. The previously unidentified term in the 
3-loop rho-parameter of the standard model is merely D$ = 6£(3) — 6Cl2(7r/3) — ^vr 4 . 
The remarkable simplicity of these results stems from two shuffle algebras: one for nested 
sums; the other for iterated integrals. Each diagram evaluates to 10 000 digits in seconds, 
because the primitive words are transformable to exponentially convergent single sums, 
as recently shown for £(3) and C(5), familiar in QCD. Those are SC*(2) constants, whose 
base of super-fast computation is 2. Mass involves the novel base-3 set SC*(3). All 10 
diagrams reduce to SC*(3)USC*(2) constants and their products. Only the 6-mass case 
entails both bases. 
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1 Introduction 



This work concerns a remarkable fusion of number, topology, algebra, and computation, 
revealed by new results for the single-scale massive 3-loop scalar vacuum diagram 

V( ri ...r 6 ):= J [dp] j[dk) j[dl) P x {k)P 2 {p + k)P^k-l)PS)P^P + l)P & {p) (1) 

where Pj(p) := l/(p 2 + m 2 rj) is a massive or massless propagator, with r| = rj, and the 
integrals are over D := 4 — 2e euclidean dimensions, with norm p 2 and a measure |J 

d D l 

[M] := m D -H D ^Y{l+e) (2) 

that removes the irrelevant constants log(47r/m) and 7 = — r'(l) as D — > 4. The symmetry 
group of (HD is that of the tetrahedron, S , 4 , whose generators are illustrated in Fig. 1. There 
are 10 distinct colourings of the tetrahedron by mass, illustrated in Fig. 2. The massive 
lines in V2A and V2N are adjacent and non-adjacent, respectively; in the dual cases, V^a 
and Van ^ it is the massless lines that are adjacent and non-adjacent; in cases V^t, V35 and 
V3L, the massive lines form a triangle, star and line, and hence the massless lines form a 
star, triangle and line. 



Fig 1: Symmetries of the tetrahedron 




Fig 2: Colourings of the tetrahedron by mass (denoted by a blob) 




v 3L V 4A V m V, v 6 

In section 2 we review previous work, which determined [l], [§ only 4 of the 10 finite 
parts analytically. Sections 3 and 4 give our method and results for all 10 diagrams. 
From them one may obtain results for all three-loop diagrams that are dominated by a 
single large mass, such as the top-quark mass 0, [5| , when the external momenta are small 
compared with that mass. Sections 5 and 6 give methods and results for the underlying 
mathematical structure || [7[]. Section 7 highlights the most remarkable of the many 
findings. 
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2 Previous results and new clues 



In each case, j, we shall evaluate and analyze the finite part of 

Vj = (i + l) 6C(3) + 3C(4) - Fj + 0(e) (3) 

obtaining simple results for the 6 cases that were previously unknown. In (|3]), 

F S := Mm(V! - Vj) > (4) 

is a scheme- independent constant, with Fx = by definition and Fj > otherwise, since 
the effect of more than one mass in Vj is to decrease the value of the positive integrand, 
in euclidean momentum space. The remaining terms in ([3D were determined from the 
complete result for V\ in D dimensions, which involves only V functions. 

Precisely three of the 10 vacuum diagrams may be reduced to V functions, using 
integration by parts for massive [|lj and massless subgraphs. A result for V± follows 
immediately from the massless two-point function of M. We gave the method for V 3 t, 
with three masses forming a triangle, in IJ. In an impressive analysis || of integration by 
parts for all 10 cases, Leo Avdeev, sadly now deceased, identified V 2 a, with two adjacent 
masses, as a third T-reducible case. The results of the present work show that there are 
no more. Performing the algebra, we obtained 



3e 4 (l 




- 2e)Vx 


6e 4 (l 




-2e)V 2A 


4e 4 (l 




- 2e)V 3T 



G3 -3G-1-1 G3-1G2, 



Gx,-2 Gx,x 
G%-x{G2-x — 3(j2,i + 2^2,2) 
Gx,x 

G3 -1G2 2 



(5) 
(6) 



Gi, 



1 (7) 



with 



G/i,u 



r(l + ne)F(l + ue) 



T(l + fxe + ue) 

= l- / iz/C(2)e 2 + ^(/i + z/)C(3)e 3 -/iz/( /U 2 + i / i I / + z/ 2 )C(4) £ 4 + 0( £ 5 ) (8) 

yielding only ^(4) = ^n 4 at 0(e 4 ). The resultant values 

Fx = (9) 
F 2A = 8C(4) (10) 
F 3T = 12C(4) (11) 

provide useful tests of the methods that we use later. 

In addition to V 3 t, we obtained in || a reduction to 3 F 2 series of V^, in the QED 
case with 4 massive (electron) lines and two massless (photon) lines that are not adjacent. 
The complete result, in D := 4 — 2e dimensions, is 

3Vw + W 3 T _ W(l, 1; 1, 0) - W(l, 0; 1, 1) 

2 e 2 (l-e){l-2e) 1 ' 
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in terms of Saalschiitzian hypergeometric series of the form [[l 

w ^ V] a ' # := — (tr^KiT^) (13) 

constrained by wreath-product symmetry. They generate polylogarithmic ladders that 
enable us to compute the ten millionth hexadecimal digits of the QCD constants £(3) and 
C(5). The expansion of flT|) yields Li 4 (|) := En>o(f)7™ 4 at 0(e 2 ). At 0(e n ~ 2 ) one 
encounters alternating Euler sums @ of weight n, which cannot be reduced to Li„(|) 
for n > 6. The result at weight n = 4 is [|I|] 

£ 4 := |(T4iv - ^st) = 16 Li 4 (i) + | log 4 2 - |tt 2 log 2 2 - ^tt 4 + 0(e) (14) 

for the finite combination that enters the 3-loop QCD radiative corrections to the elec- 
troweak rho-parameter ||, These also entail the two-loop constant defined in || as 

- g (3B+ ^ + 2) , = jj^W) = (15) 

which enters the finite part of the dimensionally regularized rising-sun two-loop diagram 
with 3 equal masses |J and involves the maximum value, Cl2(vr/3), of the weight-2 case of 
Clausen's polylogarithm C\2 m (9) := S Li 2m (exp(i6')) = J2 n >o sm(n9)/n 2m . An analysis of 
the rising-sun two-point function in D dimensions was given with Jochem Fleischer and 



Oleg Tarasov in [ID|. At zero external momentum, we obtained 



{i - a - *) ^(M; I; 1) - v Jl (e J = s> + ow He) 

with a hypergeometric series and an Euler Beta function combining to yield (|1|) as D — > 4, 
in agreement with Z)-dimensional analysis of the 2- loop vacuum diagram ||11||. 



In the course of the current investigations, we proved that S2 is, most remarkably, an 
exponentially convergent sum: 

s - 4 W 1 Vf 9 9 12 3 1 1 1 

2 81^o^ 27 J \(6n + l) 2 (Qn + 2) 2 (6n + 3) 2 (6n + 4) 2 (6n + 5) 2 J 

(17) 

which yielded 10 000 decimal places in 35 seconds on a 333 MHz workstation. 



Following [T2j we say that if the <ith digit in base b of a constant C is computable 
in time = 0(d\og T d) and space = 0(\og s d), then C E SC*(6). Then (0) implies that 
S2 ^ SC*(3), since the inverse powers of 3 are trivial in base 6 = 3 and there is an SC* 
algorithm ]12[ for the rest of the summand, with logarithmic exponents that are merely 
T = 3, for time, and S = 1, for memory. 

Mathematical interest arose from the base-2 discovery |T2[ of David Bailey, Peter 
Borwein and Simon Plouffe that {tt, tt 2 , log 2, log 2 2} C SC*(2). Very recently 0, we 
added 14 significant constants to their list. It is now proven Q that 

{7f j log k 2 I j + k < 3} U {ii 2j log k 2 I 2j + k < 5} U {G, C(3), C(5)} C SC*(2) (18) 
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where G := SLi 2 (i) is Catalan's constant, while ((3) and ((5) are intensely familiar, 
from massless 3-loop QCD. Amusingly (but not very usefully) the pentalogarithms of |J 
enabled us to compute the 64 hexadecimal digits of £(5) that begin at the 10 7 th place, 
while less than 10 4 decimal digits had been tabulated. 

It came, however, as an utter surprise that a Feynman diagram with 3 masses should 
choose b = 3 as its base of super-fast computability. In this paper, the focus is not on 
digits, in any base, but on the novel connection between the topology of the diagram and 
the type of constant in its value. We take very seriously the fact that 3 masses generate 
a two-loop constant S 2 G SC*(3), even though we have no interest in its ternary digits. 

It follows from ( |TT] , |T4| , |T8D that the three-loop diagram V^v, with 4 massive lines and two 
non-adjacent massless lines, has a finite part e SC*(2). These remarkable connections 
between the base of fastest computability and the number of massive lines will be pursued, 
vigorously, at three loops. As an example, we cite the case of the constant defined, in the 
MS scheme, by Leo Avdeev, Jochem Fleischer, Sergei Mikhailov and Oleg Tarasov as [§] 

D 3 := V 3L := lim (v 3L - (19) 

which is the only constant in the three-loop QCD corrections to the electro-weak rho- 
parameter of the standard model that had not been identified analytically, though a 
22-digit value was given in with the first 6 digits confirmed by ||. As one of many 
new results in this paper we shall show, in section 4.2.2, that it reduces to zeta values, 
together with the square of the SC*(3) constant (|P7D , as befits the 3 massive lines of V 3 l- 

It was the occurrence of Cl2(7r/3) := S Li 2 (exp(i7r/3)) at two loops that also alerted 
us to the possibility of reducing the finite parts of all three-loop single-scale vacuum 
diagrams to a (very) few weight-4 primitives of a pair of shuffle algebras involving the 
sixth root of unity A := exp(i7r/3) = (1 + i\f?>)/2. In the case of constants linked || to 
SC*(2), namely Euler-Zagier sums |7|, [13|], the operation of this pair of algebras is fairly 
well understood [|TJ] . Thanks to advice from Pierre Deligne, I was led to suppose that an 
extension to complex-valued sums involving A might result in a simple set of results for 
three-loop diagrams previously supposed too complicated for exact analysis. Thanks to 



encouragement from Dirk Kreimer, to take most seriously the idea of a connection [|15 
between the topology of a diagram and its value, I was motivated to calculate all 10 
diagrams, and found that all are reducible to primitives of the structure that is so tightly 
constrained by the pair of shuffle algebras. Thanks to David Bailey's integer-relation 
finder pslq |I6]| , I was able to relate these primitives to constants in SC*(3) U SC*(2). 



Sections 5 and 6 contain some of the large number of algebraic, analytical, and numerical 
discoveries that resulted. But first, and foremost, sections 3 and 4 give details of the field- 
theory calculations^ that both motivated and enabled the more abstract mathematics. 

In contrast to the findings, the method of calculation was prosaic [[17], [18]], requiring 
only the Cutkosky rules, Cauchy's theorem, and a large amount of perspiration. 



headers concerned with mathematical structure may wish to skip to the results in section 4.6. Yet 
to do so would be to miss out on what is involved in the discipline of perturbative quantum field theory. 
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3 Dispersive method 



Our method is strictly 4-dimensional, exploiting the dispersive results of |§. Defining the 
finite two-point function 

I(n...r 5 ;p 2 /m 2 ):=^ J d 4 k J dH P 1 (k)P 2 (p + k)P 3 (k - l)P i (l)P 5 (p + I) (20) 
in 4 dimensions, we obtain 

V( n . . . r 5 , 0) - V(n . . . r 5 , 1) = [°° dx J(n . . . r 5 ; x) {- - + 0(e) (21) 

JO IX X + 1 J 

for the difference of vacuum diagrams with a massless and massive sixth propagator. 
Suppressing the parameters r\ . . . r 5 , temporarily, we exploit the dispersion relation 



OO 



*0 



1 1 



I(x) = / dsa(s) y (22 



S + X s 



where — 2nia(s) = I(—s + iO) — I(—s — iO) is the discontinuity across the cut [— oo, — s ] 
on the negative axis. Integration by parts then gives 

/(*)= I™ ds a' (s) \ -log + log + (23) 



s L V s J \ Sq 



where the constant term in the logarithmic weight function may be dropped if <j(sq) = 0, 
as occurs when so = 0. As x — > oo, we obtain the universal asymptotic value 

6C(3) = /(oo) = / ds a\s) {log(s) - log( So )} (24) 

J sn 



with the log(so) term dropped when so = 0. The finite difference in ( pip is obtained 
from (|23|) as 



dxI{x){--^—)= rdsa'(s){L 2 (s)-L 2 (s )} (25) 
Jo Lx x + lj J so 

with a dilogarithmic weight function 

L 2 (s) := r log (^f) = Li 2 (l - 1/,) = -| log 2 ( S ) - Li 2 (l - s) (26) 

Jo x(x + l) \l+x/sj 

that is chosen to satisfy L 2 (l) = 1, thus enabling one to drop L 2 (s ) for s = and s = 1, 
which covers all the cases with N < 3 massive particles in the two-point function, and 
hence iV + 1 < 4 massive particles in vacuum diagrams. 

We now prove that the two terms in the weight function ( |26"D can be separated to yield 
the finite parts of the vacuum diagrams combined in (0), as follows: 

F(n...r B ,0) = \ cW(r 1 ...r 5 ; S ){log 2 (s)-log 2 ( So )} (27) 
F(n...r 5 ,l) = -/ dsa'{r l ...r b -s) {Li 2 (l - s) - Li 2 (l - s )} (28) 



with constant terms in the weight functions that are inert when Sq = and when s = 1- 
The proof uses the representation 

rod 

I(x) = 6C(3)+ dsa'(s){-\og(x + s) + \og(x + s )} (29) 

J Sn 



in which the asymptotic value (|24| ) is subtracted. Then one obtains 

r dx /(00) ~^ X) = - f°°ds a'(s) {Li 2 (l - s) - Li 2 (l - s )} (30) 

JO X + 1 J s 

which establishes the validity of (0), up to a possible universal constant. However this 
constant must vanish, since in the case rj = the left-hand side of (p8 ) vanishes, by 



virtue of the definition (g), and the left-hand side of ( |30| ) vanishes, because the massless 
two-point function coincides with the asymptotic value, 6£(3). The first result, (^7|), then 
follows from (^), using (p5|j26|) . It implies a sum rule for the two-point functions with a 
single mass, which we shall verify by explicit calculation. 

It is now clear that a dilog is needed to extract a finite part from a', when the sixth line 
is massive, whereas the square of a log performs the job when the sixth line is massless. 
For a vacuum diagram with N + 1 < 6 massive lines, one has a choice of method: to 
include only N massive lines in the two-point function, thus avoiding elliptic integrals || 
in a'; or to include all iV + 1 masses, thereby avoiding a dilog. We found it better to 
avoid elliptic integrals, since they are tedious to program, though speedy to compute by 
the process of arithmetic-geometric mean M. 



4 Calculation of diagrams 

There are 13 distinct cases of the two-point function (|20f) with at least one massive particle. 
We denote the spectral densities by 



(Ji(s) 


:= 


;i,o, 





0. 


0: 


5) 


03 («) 


:= 


;o,o, 


1 


0, 


0; 


5') 


0"l2(s) 


:= 


:m, 





0. 


0; 


5') 


0-13(s) 


:= 


;i,o, 


1 


0. 


0: 


s) 


014 (s) 


:= 


;i,o, 





1, 


0: 


5') 




:= 


;i,o, 





0. 


1; 


5') 


012(5) 


:= 


;o,o, 


1 


1, 


1; 


s) 


013(5) 


:= 


;o,i, 





1, 


1; 


5') 


014(5) 


:= 


;o,i, 


1 


0. 


1; 


5') 


015(5) 


:= 


;o,i, 


1 


1, 


0: 


s) 


01 (s) 


:= 


;o,i, 


1 


1, 


1; 


s) 


0^3(5) 


:= 







1, 


1; 


5) 


a(s) 


:= 




1 


1. 


1; 


5) 



where the subscripts of mass ivc indicate the massive lines, while those of maS siess show 
the massless lines. Thus (P7|,ES|) provide 26 ways to evaluate 10 finite parts of vacuum 
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diagrams, with 16 checks available. Table 1 decrypts the 26 integrals to the 10 tetrahedral 
topologies. Since only a' 15 , a[ and a 1 involve an elliptic integral, there is a systematic 
polylog route to all finite parts, save that of the totally massive case, V 6 . 

Table 1: Integral to diagram dictionary 



density 


0"1 


0"3 


<y\2 


0"13 




Cl5 


0"12 


0"13 




0"15 


0~1 


0-3 


a 


integral fl27|) 


Vx 


Vy 


V 2A 


v 2A 


v 2A 


V 2 N 




v 3L 


v 3S 


v 3L 


v AA 




v 5 


integral (p8|) 


v 2A 


V 2 n 


V 3S 


v 3L 


V 3 t 


v 3L 


v 4A 


v iA 






v 5 


v 5 


V 6 



Specializing the analysis of || to cases with r 2 = rj, we obtain 



a'in ...r 5 ; s) 



-- {^(r 1 ...r 5 ;s)e( S -(r 1 +r 2 ) 2 ) + (1^4 ) 2^5)} 
+ {a' b (n . . . r 5 ; s) 6 (s - (r 2 + r 3 + r 4 ) 2 ) + (1 <-> 2, 4 «-> 5)} 

aV ' ; i(r 4 +r 5 )2 A(s,r 1 ,r 2 ) &c Vx-s + ^O ; V ; 

^■■■^ S ) ■= m J {r3+r4)2 dx dS[ x -r 1 + i0 ) (32) 

T(s,a,b,c,d,e) := arctanh ( - _ — - _ ^ _ J (33) 

A(a,6,c) := Va 2 + b 2 + c 2 - 2ab - 2bc - 2ca (34) 

with integration by parts in ([H]) giving a logarithmic result, in all cases, and differenti- 
ation in fl32~D giving a logarithmic result when r\r 3 r^ = r 2 r 3 r^ = 0, i.e. when there is no 
intermediate state with 3 massive particles. 



4.1 Euler sums in vacuum diagrams with two massive lines 

With two massive particles in the vacuum diagram, we need only one in the two-point 
function, leading to very simple results for the derivative of the spectral density. In 
particular, the result 

sa'(s)-l f 01 "^ ' 1 ] m 

sa iy s ) - \ -n(l/ 8 ) for s G [l,oo] (6b) 

was obtained in ||, in terms of a logarithmic function 

x 

/i(x) := log(l — x) H log(x) (36) 

1 — x 

that will re-appear in other cases. It is equally simple to obtain 

s°' 3 {s) = \ $\ f TSE[ y\ (37) 
3W [ -u(l/s) for sE [l,oo] v ' 

with 

2x 

^(x) := Y^^og(x) (38) 
whose denominator will be shown to produce alternating Euler sums in F 2 ^. 
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The change of sign in (p5j 37), under the conformal transformation s — > 1/s, guarantees 
the sum rules 

roc 

= / dsa'Js) (39) 
Jo 

rod 

= / dsa' k (s) log 2 (s) (40) 

roc c _|_ 1 

= / ^(s)^-log( S ) (41) 

JO S — 1 

for = 1 and k = 3. The first confirms that <Tfc(0) = 0, allowing constant terms to be 
omitted from the weight functions of ( f24| , p7lJ?8l) ■ The second gives F\ = in fl27|) , as 
required by the definition (|J) . The third is the N = 1 case of the general sum rule M 

0= rdsa'(s){ S -±^\og(s)- S -^\og(s )\ (42) 

J s I S — 1 S — 1 J 

obtained by considering finite diagrams in which one massive line is doubled, by differ- 
entiation w.r.t. m 2 . It applies when the two-point function contains N massive lines and 
the N + 1 massive lines of the corresponding vacuum diagram, with a massive sixth line, 
are all equivalent, as is the case for V 2 a, V 2 n, V35, V 3 t, and V G . Hence the final line of 
Table 1 shows that ([42]) applies to a Xl a 3 , a 12 , cr 14 , <xi 5 and a, with the last two involving 
elliptic integrals. 

The sum rules 

6 C (3) = 2 f 1 dx lQg(x) „{x) = 2 f 1 dX lQg(x) v{x) (43) 



X Jo X 

follow from (^4|) and are easily verified analytically Less trivial are the quadrilogarithms 

- 1 dx 



r dx 

F 2A = / — {Li 2 (l-l/x)-Li 2 (l-z)} n(x) (44) 

JO X 

F 2N = / — {Li 2 (l - l/x) - Li 2 (l - x)} v{x) (45) 

JO X 



that result from (|28 



The first may be evaluated analytically by the following systematic method. Using 

- \ \oi(x) - Li 2 (l - l/x) = Li 2 (l -x) = C(2) - U 2 (x) - log(l - x) \og(x) (46) 

one may transform the integrand to products of log(x), Li (x) := x/(l — x), hi±(x) : = 
— log(l — x) and Li 2 (x) := ^(dy/y) Lii(y). Then the expansion Li fc (x) = J2 n >o xn / nk 
makes the integration trivial and produces non-alternating Euler sums of weight 4, all of 
which are proven Jl3[ to evaluate to rational multiples of 7r 4 , with 

= C(4) = C(2, 1, 1) = 4C(3, 1) = |C(2, 2) (47) 

where we use the notation of [|13j for the non-alternating Euler sum 

E f[4 (48) 

rij>n J+1 >0 j=X n j 
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which is also referred to as a multiple zeta value |T{| of depth k and weight J2j s j- 
From fl4"3j4"7|) we obtain 

F 2A = 2J q — {ilog 2 (x) + C(2)-Li 2 (x) + Li 1 (x)log(x)}{Li 1 (x)-Li (x)log(x)} 
= 8C(4) (49) 




in agreement with (10). 

We are also able systematically to evaluate all Euler sums produced by 

F 2N = 4 — {ilog 2 (x)+C(2)-Li 2 (x) + Li 1 (x)log(x)}Li (-x)log(x) (50) 

JO X 1 > 

for the finite part of the vacuum diagram with two massive line that are non-adjacent. 
Now, however, we encounter the wider R world of sums of the form 

k n n i 

= e n\ (5i) 

nj >n J+1 >0 j=l U j 

which involves a string of signs, Oj = ±1, written below the string of exponents. Such 
signs are generated by the presence of Li (— x) := —x/(l + x) — J2n>o(~ x ) n i n dH) an d 
hence (|50"D . By way of example, 

/ 1, 1, 2\ = y (-1)^ 1 (-1)"3 

= A C(4) _| C( 3 )log 2 + I C (2)log 2 2 (52) 

gives the reduction of a weight-4 depth-3 sum to sums of lesser depth, and their products. 

In 0, we discovered a simple rule for the number of basis terms for reducing all 
Euler sums of weight n: it is given by the Fibonacci number F n+ i = F n + F n _i, with 
Fy — F-2 — 1. Subsequent work with David Bailey |20| has confirmed this up to weight 
n — 11, where there are F\ 2 = 144 basis terms. Using the integer- relation search routine 
PSLQ PS , we reduced all 1024 convergent alternating sums of weight 11 and depth 11 to 
a 144-dimensional basis, finding integer coefficients of up to 30 digits, for which 5000-digit 
precision was necessary. The chance of mistaken reduction was always less than 10~ 200 . 

In the present case, at weight n = 4, there are only F 5 = 5 basis terms, which may 
be taken as 7r 4 , 7r 2 log 2 2, log 4 2, £(3) log 2 and a single irreducible alternating double 
sum. Previous experience |], with quantum field theory shows that it is particularly 
convenient to take this fifth term as 

^■-^J-^-A-l -i) (53) 

since both the contribution ( |14|) to the 3-loop electroweak rho-parameter 0, [|] in QCD, 
and also the weight-4 contribution to the 3-loop electron anomalous magnetic moment [^1|] 
in QED, are rational combinations of U^i with C(4), free of the other three basis terms. 
Thus we had a very strong expectation that the same simplification would occur in the 
present case. Indeed it does; using the methods of 0, [13j we evaluated fl50"D as 



27V 



C(4) + 8U 3A (54) 
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which is, to our knowledge, a new result. Converting this to a high-precision numerical 
value is a simple matter, since [|7| 

U 3A = |C(4) + |C(2) log 2 2 - i log 4 2 - 2 Li 4 (|) (55) 

has all || of its terms in SC*(2). 



4.2 Vacuum diagrams with three massive lines 

In the first instance, we shall evaluate vacuum diagrams with three masses as integrals 
of a dilog times a logarithm from the derivative of the spectral density of a two-point 
function with only two masses. For the massive triangle in V^t, this integral yields a 
multiple of 7r 4 . We shall show how to remove the dilog, to get a product of 3 logs, in 
the case of interest in the standard model 0, [5|], namely V^l, where the 3 masses form a 
line. The result is reducible to 7r 4 and 02(71-/3). The star case, V35, is found to be an 
intriguing combination of the other two. 

4.2.1 Zeta values when the masses form a triangle 

We begin with cr^ 4 , which was evaluated, for the purposes of QCD, in [§], with the very 
simple result 

s a ' u {s) = -2/i(l/ s) 6(s - 1) . (56) 

This affords a good test of the method, since the spectral density does not vanish at 
threshold, s = 1. One readily shows that 

4C(2) = / dsa[ A (s) (57) 



CO 



6C(3) = jf dsa' u (s)log(s) (5* 

■ ■ '- :s + 2 



(59) 

F 2A = 8C(4) = 1/ dsa' 14 (s) log 2 (s) (60) 



/oo 
rf S ^ 4 (s)Li 2 (l- S ) (61) 

with the first sum rule giving the threshold value — ai^il), and the second the asymptotic 



value 7 14 (oo). The third is the limit of sum rule (|42|) at iV = 2 and So = 1, with a constant 
term in the weight function that must be retained, on account of (0). The fourth gives 
F2A1 i n agreement with (f49D , and the fifth F 3 ^, in agreement with (|ll|). Hence we evaluate 
no new result, but rather submit the methodology to a thorough work-out. 

4.2.2 The square of Ci2(7r/3) when the masses form a line 

To obtain F 3 l, entailed by the 3-loop rho-parameter [|, |J, we calculated 

15V ; l + |s-l| V ; (s-l)(s-2) V ; V ; 
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which passes the sum-rule tests 

POO 

= / dsaUs) (63) 
Jo 

6C(3) = / dsa[ 5 (s)\og(s) (64) 
Jo 

fC(4)+8[/ 3)1 = 1/ G^ 5 (s)log 2 (s) (65) 

*/ 



2V 



with the third giving agreement with (54). There is no analogue of (|59D , since the massive 
lines in V 3 l lack the necessary symmetry. The finite part of Vzl is obtained by evaluating 

poo 

F 3L = - dsa' 15 {s)U 2 {l-s) (66) 
Jo 

which readily gives a 30-digit result for the MS finite part 

J D 3 :=6C(3) + 3C(4)-F 3L (67) 

that validates the 22-digit result in ||]. 

To identify the analytical nature of F 3 l, we remove the dilog from (|66D. The method, 
in outline, is as follows. 



1. The difficulty in (|66"D resides only in the contribution to (|62[) with branchpoint at 
s = 4; the remaining part yields only alternating sums. 

2. To separate the terms in a' l5 , one should work with the combination F 3 l — F 2 n — 
J °° dx Ixsix) / x(x + 1), where the convergence is sufficient to consider each interme- 
diate state separately, without need of preserving high-energy cancellations. 

3. By Cauchy's theorem, one may write down, by inspection, a function F 15b (x) whose 
discontinuity reproduces the contribution to a' 15 from the intermediate state in (|62|) 
with branchpoint at s = 4. 

4. Integration by parts gives this contribution to F 3 l as / °° dxl[ 5b (x) log(l + 1/x). 

5. The remaining contributions evaluate as alternating Euler sums. These contain 
C(4), and an unwonted 7r 2 log 2 2 term. 

6. Further simplification is achieved by separating from F l5b (x) the terms in 7r 2 / (x + 1) 
and 7r 2 /(x + 2) that remove singularities of its hyperbolic part at x — — 1 and 
x = — 2. It is not necessary to include these in the integral over the euclidean 
region, x G [0, 00]. 

7. The tt 2 /(x + 2) term in F 15b (^x) yields a 7r 2 log 2 2 contribution to F^ that precisely 
cancels those from the other intermediate states. 



The result is 



7Qwx TT f°° dx(x + A) l9 /a; + 2\ fx + l\ 



(x + l)(s + 2) 
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with a factor (x + 4)/(x + l)(x + 2) obtained by inspecting (|62|) at s 



The simplest hypothesis for (^g) is reducibility to weight-4 terms formed from prim- 
itives already encountered in vacuum diagrams, and their products, namely to 7r 4 , U^,i 
and S 2 = 2 4 3 _5 Cl2(7r/3). We discount tt 2 S2, since measure (H) suppresses ((2) at two 
loops, but not C(4) at three loops. To see whether the 3 candidates suffice, we define the 
5 quadrilogarithms 

dx , 9 (x + 2\ , fx + 1 



Q( n ) := / arccoslT log (69) 

Jo x + n \ 2 J \ x J 

for n G {0, 1, 2, 3, 4}, which are all values of — x at which the hyperbolic function gives a 
rational multiple of 7r 2 . Evaluating these 5 integrals to high precision and performing a 
PSLQ search, one finds that there is an integer relation between them: 

= 15 Q(0) + 144 Q(l) - 448 Q(2) + 126 Q(3) + 168 Q(4) . (70) 

Moreover PSLQ confirms that there is no integer relation, with coefficients of less than 20 
digits, between tt 2 S2 and {Q(n) \ n — . . . 3}. Since C\1(tc/3), ((A) and U^,i provide only 
3 basis terms, we must adjoin a fourth that is not reducible to these. The transformation 
x = (1 — y) 2 /y maps ( j69| ) to 

Q(n) = I' dy \og\y) log ( 1 ~ V+ £ ) # log f - ■ , y „ ■ 2 ) (71) 



o V (!-Z/) 2 J dy \l + {n-2)y + y 

with n G {0, 1, 2, 3, 4} producing only logs of y and yu — y, where y\ = 1. Noting the role 
of the sixth root of unity, A := exp(i7r/3), we chose the fourth basis term as 

Vv:= urcospwa) ^^ a, m (72) 

m>n>0 m " \ ' / 



where now the notation fl5T|) is generalized to include cases with cr| = 1, as opposed to 
merely <jj = 1 for multiple zeta values, or a 2 = 1 for Euler sums. In section 6.2.1, we 
show how to evaluate V^i to 1000 digits in a few seconds. Using 60 digits, PSLQ found 
the reductions 

Q(0) = AC\ 2 2 (n/3) (73) 
Q(l) = |Cl^(7r/3) + |C(4) (74) 
Q(2) = -Cl 2 (vr/3) + f|C(4) + |f/3,i (75) 
0(3) = -fCl^/3) + f ((4) - f t/ 3 ,i + f ^3,i (76) 
Q(4) = ±f C(4) + 8C/ 3 ,i - 8^3,1 (77) 

which were then checked to 120 digits, making the chance of misidentification less than 
10 -60 . The formulas for Q(l) and Q(2) provide an evaluation of the integral in (68) that 



is free of V31 and cancels the Us t i term from the easily obtained Euler sums, leaving 



(Ch(n/3)Y+ k(2) =m-lD 3 (7 



as the wonderfully simple sum of squares that determines the rho-parameter constant D 3 , 
previously unknown |4], |5[] analytically. 
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4.2.3 A star-triangle-line integer relation 

To evaluate the massive-star finite part 

/•oo 

F 3S = - dsa' 12 (s)U 2 {l-s) (79) 



we calculated 

f a' 12a (s) for a G [0, 1] 
<r'M = \ °m(s) for s 6 [1,4] (80) 

[ °12c( S ) for S £ [4, OO] 

with the branchpoints at s — 0, 1,4 adding successive terms. The results are 

, . , 2 arccos(l — s/2) , , 

*13«00 = r= / ^ (81) 

/s(4 - s) 



arccos(g/2- 1) log(s) 
\A(4 - a) 

/ / \ 21og(g) - arccosh(g/2 - 1) log(s) 



's(s-A) 

which were stringently tested by the 4 sum rules 

/•OO 

= / dsa' 12 (s) (84) 



o 



/•oo 

6C(3) = / dsa' 12 ( S )\og(s) (85) 



o 



oo 



= / dsa' 12 (s) S -^\og(s) (86) 



F 2A = 8C(4) = \ \ dsa' 12 (s) log 2 (s) (87) 
where the third is in agreement with (|42[) and the fourth with (^9|). 



We gave (|7"9"HS3|) to maple, which returned a 30-digit result for F 3 $. Passing this 
to the PSLQ lattice algorithm, we searched for an integer relation and were rewarded by 
another satisfyingly simple result: 

^35 = lF 3T + |F 3L = 4Cl 2 (7r/3) + f C(4) (88) 

which provides a direct relation between diagrams: 

3V 3S = V 3T + 2V 3L + 0(e) (89) 

similar to that provided by (|9l, |10| , |il"l) , namely 

3^2A = Vl + 2V 3T + 0(e) (90) 

which has pentalogarithmic corrections at 0(e). I feel that the simplicity of this third 
new result, ([89]) , 3 stars = triangle + 2 lines, is trying to tell us something important 
about the mapping from diagrams to polylogs that is provided by quantum field theory 
in 4 spacetime dimensions. It is remarkable, to put it mildly, that integration of the 
physical intermediate-state contributions (^l|,^,^), with a dilogarithmic weight function, 



results in no more than the simple mnemonic (|89|), despite the different branchpoints and 
analytical properties of the discontinuities in the various channels of the three diagrams. 
Colleagues are warmly invited to demystify the integers in (8"5|). 
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4.3 Vacuum diagrams with four massive lines 

In 0, we obtained the result ([14]) for the diagram V^n, with 4 massive lines and 2 massless 
lines that are not adjacent. We briefly revisit this case dispersively, before turning to the 
new case V^a, where the massless lines are adjacent. 



4.3.1 Integer relation, with non-adjacent massless lines 

In H we obtained the result 

*, W = ^Lt£<2e<.-4) (9i) 

in terms of the logarithmic function ([36]) and a mapping s — (1 + t) 2 /t. Transforming 
sum rule (^7|) to an integral over t, we obtain 

FiN = _ 2 /; dt amm ^ (M) _ v 4 } (92) 

The integrand can be rewritten in terms of log(t) and Lifc(±i), with k = 0, 1, 2, to obtain 
a systematic reduction to alternating Euler sums, giving 

F 4N = 17C(4) + 161/3,1 (93) 

in agreement with the QED analysis of M , where we evaluated (|T3) . It was the simplicity 



of (0), in comparison to its disguise ([14]), that suggested the importance of U 3t i in quan- 
tum field theory @ . This was confirmed by the result for the weight-4 terms in the 3-loop 
electron anomalous magnetic moment [pT]| , which are proportional ^ to 39^(4) +400L/3 i i. 
The new Euler-sum result of this paper, F 2 n = yC(4) + 8U3A, strengthens the case that 
U3A1 rather than Li4(|), is the new number to expect from Feynman diagrams entailing 
alternating Euler sums. It also provides another simple integer relation: 



AN 



2F 2N - 2C(4) (94) 



though we did not need pslq to detect this, having instead a systematic procedure for 
dealing with Euler sums. It is notable that doubling the number of massive lines also 
doubles the contribution of the irreducible alternating Euler sum. 



4.3.2 Integer relation, with adjacent massless lines 

In H, we obtained the result 



_, , . _ j -(4s - s 2 )~ 1/2 arccos(l - s/2) for sE [0,4] , , 

ai2[S) ~ i (s 2 - As)- 1 ' 2 pit) for 8 e [4, 00] ' J ' 



with a mapping s = (1 + 1) ft to the logarithmic function 



it 

pit) :=— log(*)-21og(l + t) (96) 
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which enables rapid numerical evaluation of 

("OO 

F 4A = - dsa' 12 (s)U 2 {l-s). (97) 

JO 

An integer-relation search by pslq then returned 

F iA = F 3L + |C(4) = 601^/3) + ^ C(4) (98) 
which may be rewritten as a fairly simple relation between diagrams: 

3(Vl4 - V 3L ) = 2(V 3T - V 2A ) + 0(e) (99) 
akin to (|89| , |90|) . Thus there is still no role for the suspected primitive V^i in (72), since 



both cases with 4 masses enjoy remarkable integer relations to diagrams with fewer masses. 



4.4 Vacuum diagram with five massive lines 



Applying Q2"8| ) to Q9"ip we obtain 

F . = 4 jf * e«tt!!@. {u 2 (l - <i±£) _ Li2 (-3)} (100) 

which is noi reducible to the previous 8 cases. As expected, PSLQ easily found the result, 
at 30-digit precision, when given the additional constant V 3t \ to work with, returning 

F 5 = -|Cl^(7r/3) + f>C(4) + 16^,1 (101) 

from which is absent. This reproduces (|100|) at 80-digit precision, making the chance 
of mistaken identification less than 1CT 50 . 

With this fifth new result, all envisaged constants are now in play. The key question 
is: do they suffice when the final line is given a mass? 



4.5 The totally massive case 

We were able to handle the foregoing 9 cases by methods that avoided intermediate states 
with 3 massive particles. Now there is no option, since 

roc 

F 6 = -J dsa , (s)U 2 (l-s) (102) 
involves intermediate states with two and three massive particles in 

a'(s) = W a (s) B(s - 4) + ^(s) e(s - 9) . (103) 
We may, however, simplify matters by separating these contributions in 

POO 

F 6 -F 5 = ^ dsa'(s)U 2 (l-l/s) = F a + F b (104) 

F a := ^ ds^ a (s){Li 2 (l-l/s)-C(2)} (105) 

F b := / dsa'is) {Li 2 (l-l/s)-C(2)} (106) 
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where the cancellation of the cuts at high energy is not necessary to ensure convergence, 
and F 5 is given by ( |101| ). 



The two-particle cut gives a logarithm in 



a' a {s) = — ^— \ arccosh(s/2 - 1) - ^ ) (107) 



3s(s -4) 



while the three-particle cut gives the elliptic^] integral 



^ = -9 f { dX A ^ 1 ' 1 ) x + s ~ l nn c 

b{) k x-1 A(x,s,l) A 2 (x,s,l)+xs' 1 



At large s, contributions (|107| , |108|) are each 0(log(s)/s), while their sum is 0(log(s)/s 



2^ 



The integrals ( |105| , |106| ) converge separately, thanks to the ((2) in their weight functions, 



to which the combination (|104 ) is blind. 

High-precision numerical evaluation of ( [L05|) presents no problem, since it is of the 
same character as previous cases, involving merely integration of the product of a dilog 



and a log. On the other hand, it appears from (|106| , |108|) that we are also obliged to 



integrate the product of a dilog and an elliptic integral. In fact, this is not necessary, 
since we may use the method of || to circumvent such a contingency, by reversing the 
order of integration. Setting x = 1/u 2 G [4, oo] in ( |108| ), which now becomes the outer 
integration, and s = (l/u + v)(l/u + 1/v) G [(l/u+ l) 2 , oo] in the inner, we then integrate 
by parts on v G [0, 1] to convert the dilog to a product of logs, with the result 





- 2 S 




Jo 


A(u) 


■= log 


B(u,v) 


:= log 


C{u,v) 


:= log 


D(u,v) 


:= log 




' dB(u, v) 



l-u 2 
(1 + uv){u + v) 

U + V + uv 2 
(1 + uv){u + v) 

U 2 V 



dv C(u, v) D(u, v) (109) 

\ ov 



(110) 
(111) 
(112) 
(113) 



'1 + 2uv + v 2 + (1 - f 2 )Vl -4m 2> 
. 1 + 2uv + v 2 - (1 - v 2 )Vl -4m 2 , 



which establishes that i*& is the integral of a trilogarithm. As ( |1 1 1| - 4TT3| ) are logs of rational 
functions of v, the inner integral gives terms of the form Jq 1 log(c + v) \og(d + v), with 
parameters {b, c, d} that depend on u. It is shown in section 8.4.3 of |22[] that all such 
integrals evaluate to Li 3 and products of logs and dilogs. Hence, in this apparently most 
difficult case, we still arrive at the integral of a trilogarithm, as in previous cases, though 
until now the trilog factored into products of logs and dilogs. 

At this juncture of the investigation, we face a genuine quandary. If the totally massive 
case is as complicated as the Cutkosky rules suggest, then the route is clear: one should 



2 I am told that Kallen was disappointed to find that the two-loop electron propagator jl8| involves 
an elliptic integral, unlike the simpler photon propagator jL7|. 
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take the trouble to transform (|108|) to the real parts of complex complete elliptic integrals 
of the third kind and evaluate the integrand of ( |106[ ) by the process of arithmetic-geometric 
mean, as was done in ||, with a simpler mass case and a simpler weight function. If, on 
the other hand, one follows the intuition that lives in the same universe of quadrilogs 
of the sixth root of unity as the 9 previous cases, one has a lower cost (and higher risk) 
route to pursue: to evaluate the double integral ( |109| ) to modest precision, without further 
analytical effort, and then hope to discover a simple integer relation to previous cases. 
Emboldened by discussion with Dirk Kreimer, I took the latter route. 



In this optimistic strategy, the transformation achieved in ( |109[ ) is especially fortu- 



nate, since the nag routine D01FCF is notably efficient at evaluating rectangular double 
integrals in double-precision FORTRAN, i.e. to almost 15-digit precision. Adding the easily 
computed integral ( |105| ) and the already known finite part (|101| ), we thus quickly arrived 



at a numerical value of Fq that is rock solid to 14 digits. Such precision is rather modest, 
in comparison with previous cases, yet amply sufficient to discover that 

F 6 = F 3S + F 4N -F 2N = 4 (CI^tt/3) + 4((4) + 2Z7 3j i) (114) 

from which V^i is absent. This corresponds to a direct relation between diagrams 

V 6 + V 2N = V 3S + V 4N + 0(e) (115) 

verified to 15 digits. It stands as testament to the oft remarked fact that results in 
quantum field theory have a simplicity that tends to increase with the labour expended. 

4.6 Tabulation of results 

For ease of reference, we collect the 10 results in Table 2, which gives the MS finite part 



V j := lim (Vj - = 6C(3) + 3«4) - F 3 



e j 



= 6C(3) + z i C(4)+^C/ 3 ,i + SiC1^7r/3)+^y 3l i (116) 

by specifying non-zero coefficients of the 4 basic quadrilogarithms. In each case a numer- 
ical value is appended; section 6.3.3 extends it to 10 000 digits in seconds. The first and 
last values spell out the efficacy of the MS scheme: the largest absolute value is only 0.2% 
greater than the absolute minimum, \{V\ — Vq), attainable in any MS-like scheme. 
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Table 2: MS finite parts 





Zj Uj Sj Vj 




v x 


3 


10.4593111200909802869464400586922036529141 


v 2N 


—o 

-f -8 


H 0AA70C0i:A/l AH 07/17IT/IO'l OA H A/1 OCO/^OOCA/I OAVH C A 

1.800/252504018 J 4f 54818410486362860430 J 161 
1.1202483970392420822725165482242095262757 


T T 

V 3T 

V 3S 

T/ 
V 3L 


— y 

- T -b 


O COOrr , 7rO A A /IAA70A'I h O A C/^A/IOAAOAH oh h h oaoooo 

—2.52856/6844426 ( 80 1124560429980 18 11 1803828 
— 2.obUoDzz241oyoz r obUz r z ro45o r ( foz4iyi/bol4 
— o.vzi uuy'iyoyo i oozuiy i ooo/ 'ii ui / ooyoi /odioui 


v m 


_ZZ _ 6 

12 u 

-14 -16 


-5.9132047838840205304957178925354050268834 
-6.0541678585902197393693995691614487948131 




469 8 1 fi 

27 3 


-8.2168598175087380629133983386010858249695 


K 


-13 -8 -4 


-10.0352784797687891719147006851589002386503 



5 Shuffle algebras with the sixth root of unity 



It may come as a surprise that we reduced the 9 non-zero values of Fj to merely the 4 
terms {((4), U 3>1 , CI^tt/3), V 3>1 }, with 5 integer relations, given by (IggySqpjSgyil^. The 
simplicity of the results for the Feynman diagrams reflects the tight constraints that are 
placed on 

k \pjUj 

= e n 

nj>nj + i>0 j=l 





n/ 



(117) 



by a pair of shuffle algebras entailing A := exp(27r/3) = (l + iy/3)/2. In the restricted case 
Pj G {0,3}, giving X Pj = ±1, both shuffles algebras for the corresponding Euler-Zagier 
(EZ) sums have been intensively studied 0, [13], 0, [19], H HI ll> H- The 



extension from the EZ case to the somewhat harder case with pj £ {0,1,2,3,4,5} was 
rather clear, in principle []I4" |, though little seemed to be known about the enumeration 
of primitives. In sections 5.1-5.3, we consider all of the sixth roots of unity. Then we 



remark on the context BO] in which the primitive root, A, plays a special role. 



5.1 Depth-length shuffles 



Considered as a nested sum of depth k, (|117|) has an obvious property: the product of 
two nested sums, of depths k\ and k\, is a linear combination of nestings at depth k\ + k 2 , 
with some additional terms to take account of equality of summation variables. Thus the 
simplest case, with depths k\ — k 2 — 1, is 



r 1 m 2 = m 2 \+r 2 ' 1 41T * 2 niai 





with the terms on the right corresponding to the summation regions with ri\ > n 2 , rii < n 2 , 
and rii = n 2 . The final term has the same weight, but lesser depth. The generalization 
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to the product of sums of arbitrary depth is clear: there is a shuffle of the depth-length 
strings of argument pairs (sj,pj) and extra terms, of lesser depth, in which subsets of 
arguments Sj are added, with the corresponding powers Pj of A added mod 6. 



5.2 Weight-length shuffles 



The second type of ring structure results from expressing ( |117| ) as an iterated integral |3T 



In the present case, this may be regarded as a word formed by concatenation of letters from 
the 7-letter alphabet A\ := {Q, u , u)\, u 2 , uj 3 , uj^, u 5 } consisting of the 1-forms fl := dz/z 
and uip := dz/(\~ p — z), with uj p+ q = ui p . The concrete example in ( |T2"D serves to illustrate 
the equivalence of nested sums and iterated integrals: 

J 3, 1 \ f 1 dz x f*i dz 2 dz 3 f*a dz 4 2 

\ A , A J Jo Z\ Jo z 2 Jo A 6 — z 3 Jo A a — z^ 

with the third form standing as a mnemonic for the second. To prove (|119 ) one expands 

dz 



^, : - =nJ2x pr z r (120) 

r>0 



p ■ X-p-z 



which makes the Zj integrations trivial and gives 

\3r+5s X 2n 2 / o 1 \ 

c(" 2 ^ 5 ) = £ T^pi = T, W = c U?' y) < 121 > 

r,s>0 V + S ) S ni >n 2 >0 n l n 2 \ A > A J 

with rii = r + s and n 2 = s. Thus (|119| ) may be regarded as ((W) with W = Vl 2 u 3 u^ 
specifying the corresponding iterated-integral word. It is purely a matter of book-keeping 
to transform a word to a nested sum, and vice versa. The general dictionary is 

c(^-h 1 ^ 1 ^ 2 ---^HJ:=c( A S; %- P1 ■ ( 122 ) 

showing that the length of the word is the weight, n := J2j s ji °f the sum, and the number 
of Q's in the word is n — k, where k is the depth of the sum. 

The second shuffle algebra is merely the observation that 

((W 1 )((W 2 )= £ C(Wj) (123) 

where the sum is over all words Wj in the set Si j2 obtained by shuffling W\ and W 2 , while 



preserving the order of each. This is a property of all iterated integrals |JI . 
As an example at weight 4, consider 

C(1M) C(^s) = OPoJiShJs) + 2 C(fi 2 wiw 5 ) + (wi «-> cj 5 ) (124) 
which translates to a rather non-trivial double-sum identity 

£ + C>^/3) = ( + ^) (-l)*c« + l™) (125) 
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using Y,n>o cos(7m/3)/n 2 = (vr/6) 2 from p2| , on the left, and translating words to double 
sums, on the right. 

It is clear that such weight-length shuffles cannot yield the whole story, because the 
depth-length shuffle ( |II8| ) has a final term in which the number of f2's increases, which 



can never occur in ( |I23| ). Thus one should solve for both ring structures. There is a 



systematic method for solving either in isolation, using Lyndon words [26]. These form a 
subset of words, defined by a lexicographic ordering criterion, with elements and products 
that furnish a basis to which all words can be reduced by a shuffle algebra. Choosing a 
set of Lyndon words is equivalent to choosing a minimal set of brackets that solves the 
Jacobi identity in a free Lie algebra. In our case, the corresponding free Lie algebra has 
7 generators, so such a systematic technique is particularly welcome. However, solving a 
second shuffle algebra for the Lyndon words of the first appears to be a matter of brute 
force. And even then we are not done. 



5.3 Transformations of words 

There are yet further constraints, beyond those obtainable by either type of shuffle. Con- 
sider, for example, a word formed from the 3-letter EZ sub-alphabet A 3 := {Q,u ,u 3 }. 
Let us denote by ((£1; u ; u 3 ) an arbitrary sum corresponding to some word W all of 
whose letters^ are in A 3 . Its formal arguments are the three possible letters, which we 
separate by semi-colons, for emphasis. Now we transform the path from z — to z — 1, 
in the iterated integral, by the mapping 

1 — z , 

z -> (126) 

l + z K J 

which interchanges the endpoints, z = and z = 1. This induces a linear transformation 
of the 1-forms Q := dlog(z), Uq := — dlog(l — z) and uo 3 := — <ilog(l + z) from which the 
integral is constructed. Hence we immediately arrive at a highly non-trivial result: 

C(fi; W3) = COo - ^3; + ^3; w 3 ) (127) 

in which the formal arguments are transformed, and the tilde is the instruction to write 
the resultant word backwards, on account of the interchange of endpoints. By way of 
example, we apply ( |127| ) to 

U 3A := C(^Wo) = C ((w - uj 3 ) 2 uj 3 (fl + us)) = C ((O + w 3 ) ^3 (w - uj 3 ) 2 ) (128) 

= ((SkwJfi) - C(^3^W3) - ({towlwo) + C(^MD 

+ C(o&>o) - C(^ u 3 ) - ({ulu ) + ((ul) (129) 

which turns a double sum into a combination of sums of depths 3 and 4. 

There are further constraints. From the complex conjugation A — » A -1 we obtain 

C(0; ^0; u X ] u 2 ; u 3 ; u 4 ; u 5 ) = ((Q; u ; u 5 ; u 4 ; u 3 ; u 2 ; Ui) (130) 
3 With mild abuse of notation, we later abbreviate this cumbersome locution by writing C(W0 S A3. 
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where the bar denotes the complex conjugate of the sum. We can also transform z — > 1— z. 
Acting on multiple zeta values (MZVs), this gives the duality relation 

C(fi; cj ) = CK; n) (131) 

with £(4) = C(2, 1, 1) in (TQ) as a simple example. In A\ := {fi, cj , cji}, it generalizes to 

C(fl; u ] ui) = C(^o; ^; -^l) (132) 

since A — 1 = —A -1 . We find that ( |132j ) is responsible for extra reductions within the 
complex- valued sub-alphabet At, just as extension of the Zagier duality ( |131| ) to trans- 
formation (|127|) in the Euler-Zagier sub-alphabet As reduces many real-valued EZ sums. 

Next, z — ► z 2 acts on MZVs to give EZ sums: 

C(fi; uj ) =C(tt + tt; uj + lu 3 ) (133) 

which generalizes to the duplication formula 

cu ; ^ 4 ) = C(f3 + fi; c^o + ^3; ^1 + ^4; ^2 + ^5) (134) 

acting on ^4 A 2 := {Q, u , u 2) 004}, i.e. on cube roots of unity. Similarly, z — > z 3 acts on 
square roots of unity in A3, to give the triplication formula 

C(Q; u ; lu 3 ) = ((Q + Q + Q; cu + uj 2 + tu 4 ; t^i + u 3 + u 5 ) . (135) 



All known relationships between EZ sums follow from combining (|127|) with the two 
shuffle algebras and the simple relations (|13iyi33|) . Moreover there is information in ( 127|) 



that cannot be extracted from the other four sources. A case in point is the reduction 0] 

( 1 \m+n 

U^:= E M-T- = SC(6) - IC 2 (3) (136) 



m>n>0 



which is impossible to deduce from the machinery of [23], but follows by augmenting it 



with ( |127| ). It is believed that the Fibonacci enumeration of EZ sums in 0, with F n+ i 
basis terms at weight n, follows from what is given above. This is proven, by brute force, 
for weights n < 8. At n = 8, the nexus of 5 constraints reduces 3 8 = 6561 words to a 
basis of size merely F 9 = 34. Moreover, pslq reduced all 1024 finite words of length 11 
in the sub-alphabet {uj ,u 3 } to a basis of the expected size F 12 = 144. 

In all of this, one needs a regularization of those iterated integrals (less than 27% of 
the total) whose words begin with ujq, or end with Q, since these give divergences at z — 1, 
or z = 0. They are handled rather neatly by reduction to a Lyndon basis of the weight- 
length shuffles, with product terms separating off divergences and sub divergences in a 
manner analogous to the subtraction procedures of quantum field theory. The primitives 
that occur for the first time at weight n are then in the finite parts of divergent words of 



length n, again in the manner of quantum field theory |28L 32 . 
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5.4 A perspective on roots of unity 



While practical calculation entails all of the tools above, there is both physical [15[ and 



mathematical pD| interest in restricting attention to primitive terms, i.e. those terms in 
the basis at a given weight that do not factorize into terms of lesser weight. From this 
point of view, one may ignore the left-hand side of ( |123| ) and solve the weight-length 
shuffle, with great economy, in terms of Lyndon words. 

In the EZ case, one may conveniently frame |30| such a study, via the Campbell-Baker- 
Hausdorf route, in terms of a (suitably regularized) path-ordered exponential 

r 1 h ( dz dz dz \ ,->r,~\ 

7(e ,ei,e_i) = P exp / — e h e x + e_i — — (137) 

jo 2m \ z z — 1 z + 1 J 

where {eo, ex, e_i} are non-commuting objects that generate a free Lie algebra. Then the 
Lyndon words correspond to a minimal choice of Lie brackets, consistent with the Jacobi 
identity. To specify, say, the U^x content of all weight-4 EZ sums, it is sufficient to give 
the rational combination of brackets that multiplies T^U^^ in the expansion of (|137 ), since 



products of brackets are irrelevant to primitives. The question then arises: can one arrive 
at the primitive content without using the depth-length shuffles that are hidden by the 
integral representation? 

Clearly, if one is going to neglect, pro tempore, the powerful information of section 5.1, 
and has already fully exploited, via Lie brackets, the information of section 5.2, the trans- 
formations of section 5.3 are crucial. Let us denote by Xj(e , ex, e_i, e ro ) the combination 
of brackets that multiplies a primitive term, such as (ih) 3 ((3), in the expansion of ( |137| ). 
Here eoo := — eo — t\ — e_i is an auxiliary argument that serves to expose the underlying 



structure. Then ( |127| ) tells us that 

= Xj(e , ex, e-i, e^) + Xj(ex, e , e», e_i) (138) 

which involves two elements of the symmetry group, -D 4 ~ Z 2 I Z 2 , of a square. 

So far we have gained no new information, since (|138|) is merely the consequence of 



transforming the path in ( |137| ) by z — > (1 — z)/(l + z). The remarkable finding [RJJ of 
Pierre Deligne is that there is a second key property of Xf 

= Xj(e , ex, e_i, e 00 )+X j (e l , e ro , e , e^+X^e^, e_i, ei, e )+X j (e^ 1 , e , e^, e x ) (139) 

in which one sums over four elements of D^. 

We have verified ( |139| ) in the case of the brackets multiplying C(3), using the complete 
results from previous sections. Unlike ( |138| ), the new constraint ( |139| ) is not a general 
property of EZ sums, but only of those parts of them that do not involve reductions to 
products of terms of lesser weight. In thus helps to predict the coefficients of new terms, at 
a given weight, while those of the products of old terms are ignored, in the firstfj instance. 
In this sense, it sits very close to the preoccupations of K71, 123, R2L RBI R3, RR RH], where 



Bob Delbourgo, John Gracey, Andrei Kataev, Tolya Kotikov, Dirk Kreimer and I were 
concerned with the constants that first enter quantum field theory at a given order in 
perturbation theory. 



'Pierre Deligne has devised p3| an iterative Lie-algebraic procedure for progressing beyond this stage. 
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Now the crux of these observations is that when Pierre Deligne communicated ( 139 ), 
he included the following thought-promoting remark |30] : 



If A = | ( 1 + V~ 3 J (sixth root of 1), one could hope for having a similar 
story for the free Lie algebra generated by e , t\ and e\, and for 



„ f 1 1 / dz dz dz . 

7 = PeXP io 2^( v e °- + ei — +6A — 1 • (1R,) 



z z — 1 z — A 



Putting this remark together with the commonplace finding of Cl2(7r/3) := 3Li 2 (A) at 
two loops ||, we undertook the task of identifying the relation of D 3 := V$l, in the 3-loop 
rho-parameter, to polylogs of A, expecting to find a new weight-4 primitive, to account 
for previous difficulty |4], |5| in identifying this term. As seen in Table 2, the situation was 
far simpler: none occurs; we encounter only products of sums of lesser weight. Moving 
up to 4 masses, we expected to detect at least V^i. As seen in Table 2, the situation was 
equally simple, with reductions to cases with fewer masses. With 5 masses, V^i at last 
appeared; yet with 6 it was again absent. Such simplicity in quantum field theory seems 
to be telling us that the primitives are very few in number, at weights n < 4, and that the 
topology of the diagram determines which of those few appear in the quantum amplitude. 

Accordingly, we shall conclude the paper with a study of all words with less than 4 
letters, and of all 4- letter words with depth k < 2, as entailed by the Feynman diagrams. 
We shall pay particular attention to the Deligne sub-alphabet, A\ := {Q, uj$, ui}, with the 
generalized duality relation (|132| ). The physical motivation is that Table 2 implicates Ai in 
the especial simplicity of 3-mass diagrams, while the more familiar EZ sub-alphabet, A3 : = 
{Q, a;o, uj 3 }, enters with 2 masses. We show that these sub-alphabets choose different bases 
of super-fast computability: 6 = 2 for Euler-Zagier primitives; b = 3 for Deligne primitives. 



6 Poly logarithms of the sixth root of unity 

The convergent one-letter words of the 7-letter alphabet are 

i i-7r for p = 1 

C(w p ) = -log(l- A p ) = <| - 1 log 3 + %\r forp = 2 (141) 

— log 2 for p = 3 

with p —>■ 6 — p for complex conjugates. Thus 3 primitives are presumed at weight^ n < 2. 



5l 



Physicists tend to set Ti := h/2n = 1 in (137); mathematicians tend to put h = 1, to rationalize 



powers of 7r. Thus the weight of n is a matter of convention. All parties agree that it is less than 2. 
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6.1 Reduction of all two-letter words 



At the level of dilogs, one obtains the following single sums from |£2 



((Qu p ) = Li 2 (A p ) 



¥ 2 

> 2 



+ 
+ 



for p = 
zC1 2 (tt/3) forp = l 
|C1 2 (tt/3) forp = 2 

for p = 3 



(142) 



where Cl 2 (7r/3) := J2n>o sin(7m/3)/rz 2 is presumably primitive. For the double sums, 
((u p u) r ), we need consider only the 9 cases of Table 3. 

Table 3: Real and imaginary parts of 2- letter double sums in the 7- letter alphabet 
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3? ((UpUJ r ) 
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V/3) 










2 





72 " ^ 


| log 2 3 
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V/3) - 


^TTIO 
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|log 2 2 
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^Trlog 


3 
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1 




l Li 2 (i) - 


|lo { 


5 2 2 + ilog21og3 


-|ci 2 


>/3) - 


|vrlog 


2 + 


ivrlog 


;3 


4 


1 


~^ 2 + 


|Li 2 (|) + 






-|ci 2 


V/3) + 


l^rlog 


2 - 


ivrlog 


;3 


5 


1 










-|C1 2 


>/3) + 


fvrlog 


3 






3 


2 




|Li 2 (|) + 




5 2 2 


±C1 2 


>/3) - 


^Trlog 


2 






4 


2 


^ 2 + 


±log 2 3 






-±C1 2 


>/3) + 


^7Tl0J 


?3 







The remaining cases are then determined by the shuffle 

log(l - X p ) log(l - X r ) = CKo; r ) + C(u r u p ) 



(143) 



and by the complex conjugation (|130|) . Making extensive use of identities in |22|, we found 
that there is a second weight-2 primitive, which may be taken as Li 2 (|). 

Fast computation of the primitive Li 2 (^), in the real parts, is straightforward. Almost 



as fast, is the evaluation of 



Cl 2 (vr/3) = 



dylogjy) 
o 1 - y(l - y) 



^ 2n-l ^ 



i>0 n 



( 2 ;) 



k=n 



k 



(144) 



in the imaginary parts. Expanding the integrand in powers of y(l — y), we obtain a 
summand that is a derivative of an Euler Beta function. Computation of the resulting 
inner harmonic sum, in the final double sum, is performed in the same loop that handles 
the outer sum, which is exponentially convergent, having the same factor of l/4 n as 
Li 2 (|). From a computational perspective, the two primitives are of the same highly 
convergent character. It thus takes only a few seconds to compute all 36 convergent two- 
letter words to 1000-digit precision, and hence to obtain to this precision the reducible 
quadrilogarithm (|78| ) in the rho-parameter. In section 6.2.5 we shall go a step further, by 
reducing v^3Cl 2 (7r/3) to an exponentially convergent single sum. 
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It is notable that ((uj 3 uji) of Table 3 entails 7 of the 8 real basis constants at weight 
n = 2. We conclude that care should be taken when restricting attention to the 4- 
letter union A\3 := A\ U A3 = {Q, u , u)i, u 3 } of the Deligne and EZ sub- alphabets. 
Though such a restriction might appear tempting from the perspective of the weight- 
length shuffle ( |123|) , it has little to commend it from the perspective of dictionary (|122|) , 
which shows that it denies access to the full content of depth-length shuffles. 

A comparable situation occurred in M, where we showed that there are features of 



MZVs, in A := {f2,a>o}, that are better [271, |28] understood by extending the analysis to 



EZ sums, in A3 := oo 0} 003}. Notable among these was the discovery |7| at weight 12 
that the depth-4 multiple zeta value £(4, 4, 2, 2) := ( ((Q 3 lv q ) 2 (Qlj ) 2 ) G A is not reducible 
to terms of lesser depth in A, yet is reducible to the double sum Ug t s := ((Q s U3Q 2 lj ) G 
.A3. Such phenomena, which we refer to as "pushdown", make Don Zagier's challenge of 
enumeration of primitives of A, bigraded by weight and depth, a much tougher task than 
the enumeration of primitives of A3, conjectured in M. Hence the conjectured answer for 



reducing MZVs to MZVs, given with Dirk Kreimer in |28[ , is of a much subtler form than 
those for reducing EZ sums to EZ sums, or MZVs to EZ sums. 

Heeding such experience, we are prepared for reductions of sums in A13 to products of 
primitives of the full alphabet, A\. Inspection of the results above reveals that this has 
already occurred, since they imply the amusing evaluation 

sin (|vrm — ^Tfnj := 3 {((ujiuj ) — Q((ujiuj3)} = tx log3 . (145) 

m>n>0 mn 

Remaining strictly within ^4.13, we would be in the ironical situation of accounting 7rlog3 
as primitive, whereas we know from (|141|) that in A\ it reduces to — 63 {C(^i) C( w 2)}, 
reminding us that A 2 , the primitive cube root of unity, is also a sixth root. 

Notwithstanding the note of caution sounded above, there is still a lively hope that the 
Deligne-Euler-Zagier (DEZ) sub-alphabet A13 is of the essence in quantum field theory, 
since the results for all 10 vacuum diagrams of Table 2 lie within it, entailing only 

C(4) = an 3 u )eA (146) 
U 3;1 = ((flVo)ei 3 (147) 

c\i(n/s) = ^ a^)+a^o) J eAi (148) 

V 3A = »((OWi)£ii3. (149) 



6.2 Reduction of all three-letter words 

We already know, to 1000-digit precision, the subset of 16 product terms 



^prod 


= (H 2 U1 2 )U(% 3 UI 3 ) 




(150) 


n 2 


= {7rCl 2 (7r/3), 7r 2 log2, 7r 2 log3} 




(151) 


12 


= {C1 2 (tt/3) log 2, C1 2 (tt/3) log 3, vr 3 } 


(152) 


n 3 


= {Li 2 (|)log2, Li 2 (l)log3, log 3 2, 


log 2 2 log 3, log 2 log 2 3, log 3 3} 


(153) 


i 3 


= {vrLi 2 (i), 7rlog 2 2, 7rlog 2 3, 7rlog 


; 2 log 3} 


(154) 
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in the basis for reducing real and imaginary parts of the 252 convergent three-letter words 
of A\. In accord with remarks in the previous section, we also expect these products in 
reductions of the DEZ sub-alphabet, A13. However, we know |7|] that only the products 
7r 2 log2 and log 3 2 occur in the words of A3, all of which are real, and we expect only 
7r Cl2(7r/3) and 7T 3 as products in real and imaginary parts, respectively, of words in A\. 

Note that in (|150|) we have partitioned the product terms into those that are known 
to occur in real (TZ) and imaginary (X) parts, and further partitioned them into IZ2 UZ 2 , 
known to occur in reductions of three- letter double and triple sums, and TZ 3 U Z 3 , which 
can not be produced by the shuffle algebras in the case of double sums, but certainly 
occur in reductions of triple sums. 

To these products, we must adjoin the primitive C(3), in the real parts of 



C(n 2 cu p ) = Li 3 (A") = 
and their complex conjugates, with p - 

u a , b := ((n*- 1 ^- 1 



C(3) 
K(3) 



+ 



-|C(3) + 



L 162" 
i-TT 3 



-!C(3) 



for p 
for p 
for p 
for p 




1 

2 
3 



(155) 



6 — p. Turning to double sums, we define 

l ]\m+n 



^0 



E ■ 

m>n>0 

E ■ 

m>n>0 



m a n b 



cos(27rn/3) 



m a n b 



(156) 
(157) 



where U a ^ is proven to be reducible at weight a + b = 3. We shall show that the same 
applies to V a ,b at this weight. However, close inspection of Eq (8.111) of |22] reveals that 



Li 3 (z/^3) 



1 



E 

n>0 



■1/3)' 



i (-1/3)" 
V3^o(2« + l) 3 



(158) 



will be generated by some of the three-letter words, as will 



Li 3 (A/2) 



E 



n>0 k 

iy/3 



+ 



2(-l/8) n (-l/8) n (-l/8) n ' 
(3n + l) 3 ~ (3n + 2) 3 ~ (3n + 3) 3 

-1/8)™ , (-1/8)™ 1 



E 

n>0 



+ 



(3n + l) 3 (3n + 2 



(159) 



These exponentially convergent polylogs are readily computed to high precision, yet pslq 
found no relation between their real and imaginary parts and the product terms already 
assembled. Hence we adopt 



Swork := Sprod U (C(3), Li 3 (i/V3), Li 3 
working set of constants at weight n = 3. 



(160) 
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6.2.1 High-precision reduction of trilogarithms 

We can achieve rapid computation of 



^ = - 5 y Mi-y + y) = -E^ (i6i, 



2/ n>0 

,2> 



by expanding log(l — y + y ) in powers of ?/(l — y), to obtain 

1 r i 2n-l ~| 

C n = -2CU + - + V - (162) 




with Cq = C(2). Hence we arrive at a fast numerical algorithm that updates the harmonic 
sum, the central binomial coefficient, the value of C n , and the truncated value of V^i, in a 
single loop, with the same exponential convergence as for (144). There is one point of note: 



at intermediate stages one needs to handle 50% more digits than the target precision, since 
recurrence relation ( |162|) magnifies any error in the starting value, Cq = |vr 2 . Generalizing 



this method, we also obtained a 1000-digit value for V^i in a few seconds. 
At 100-digit precision, pslq easily achieved the reduction 

V 2A := & ((QoJsOJi) = -§ C(3) + js log 3 + § it C1 2 (tt/3) (163) 

which we then checked at 1000-digit precision. Indeed, all searches for reductions of the 



real parts of three- letter double-sum words produced only £(3) and products in (|151|) , 
giving no cause to enlarge the working set ( |160| ). 

Moving on to the imaginary parts of double sums, we evaluated to 100 digits the 4 
integrals with a 2 = t 2 = 1 in 

W„ T := - f 1 dyl0M arctan ( -^L) ( 164 ) 
Jo 1 + ay \ 2 — ry J 

and obtained the PSLQ reductions 

QfC(^oWi) = = 3I4 7T 3 (165) 

-9fC(fiw 3 Wi) = = -3§j7r 3 + |{QfLi3( 4 yV3)-> log 2 3} (166) 

3fC(^ow 2 ) = WCi,_ 1 = 5 ^7r 3 -ia 2 (7r/3)log3-|W 1 , 1 (167) 

-9fC(^ 3 w 2 ) = ^x^^-iC^TT^logS-W!,! (168) 



showing a simple product in (165) and an anticipated primitive in ( |166| ), with the combi- 
nation SLi 3 (i/\/3) — ^g7rlog 2 3 common to ( |166| - |168|) . Hence it is better to use ( |166| ), in 



lieu of the imaginary part of ( |158j ), for integer-relation searches. We cite 

S((^ 3 ) = 4vr 3 -|Cl 2 (7r/3)log2 + H/ lil -2{SLi 3 (A/2)-^7rlog 2 2}(169) 
9fC(fiwiw 2 ) = il^ 3 - Cl 2 (vr/3) log3-2Wi,i (170) 

as two more findings of PSLQ that are significant for later work. The first shows the 
combination 3 Li 3 (A/2) — ^7rlog 2 2 that enters the imaginary parts of some double sums. 
The second is notably free of this combination. 
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6.2.2 One primitive trilogarithm in Ai 



In the light of ( |165| ), it becomes rather interesting to see if the Deligne sub-alphabet, Ai, 
has C(3) as its sole primitive trilogarithm. According to pslq, there are no more. We 
shall prove this, by establishing the reductions^ 



C(3) 



7T 



C1 2 (tt/3) 



1 324 " 



C(3)-|7ra(7r/3) + i^7r 



2 
3 

i gl 7i 



(171) 
(172) 
(173) 



which show that these 3 sums and the remaining 7 



(174) 
(175) 
(176) 
(177) 
(178) 
(179) 
(180) 

have, at most, ((3) primitive term. 

We have reduced the problem to proving (|171| - |173| ), by using generalized duality and 
weight-length shuffles to determine the remainder. The three extra relations are 

C(^V) 











= CK)C(^i) 


- 2C(ficj 1 2 ) 


C(^i^M)) 


= CK)C(^q) 


- ((fiu^o) - C(^o^i) 


C(^i^i^i) 


= too 




C(^iu;iu;o) 


= CM 










C(^icjo^o) 


= -c(fiV) 





C(wi)C(nwo) 

C(nw cui) + C(fiu;iu;o) 
C(wiflwo) 



C(cJificJi) + £(fiu;ou;i) 
-C(^o<^i) 



(181) 
(182) 
(183) 



The first is a depth-length shuffle. The second is obtained by equating the depth-length 
and weight-length shuffles for the product ((uo)({Quji), from whose difference the diver- 
gent term then cancels. The third is a further example of generalized duality. Solving the 
10 equations ( |174j - TS3 ) we determine the 10 convergent sums of depth k > 1 in terms of 
singles sums and their products. Then Leonard Lewin's book p2| provides the proof that 
the p = 1 case of ( |155| ) involves no new primitive, leaving £(3) in the p = case as the 
sole candidate^ for a primitive trilogarithm in Ai := {Q,uj ,uji}. 



6.2.3 Algebraic reductions in A\ 

Here we give an algebraic method that determines 123 of the 126 convergent double sums 
in the real and imaginary parts of three- letter words of the full alphabet, A\. The tally 
of double sums comprises 36 + 32 from real and imaginary parts of ((Qui p u) r ), and 30 + 28 
from ((ujpflUr) with p ^ 0. We are able to reduce all 66 real parts and 57 of the 60 
imaginary parts, by systematic computer algebra, as follows. 

6 They were all found by PSLQ, long before we found their proofs. 

7 That C(3)/7r 3 is not a rational seems to be as far from the possibility of proof as it is from doubt. 
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The first step is to use the weight-length shuffles 



((u> p )((Qu) r ) = ((UpQuJ r ) + ((Qu p LJ r ) + ((Qu r UJ p ) (184) 

to eliminate all instances of ((u> p Qu> r ), including those with p = 0. Then one may use all 
36 of the complex depth-length shuffles 

(((jjp)((ttuj r ) = ((ujp£lLUp +r ) + ((Qu r LU p+r ) + ((£l 2 LUp +r ) (185) 

since in those with p = the divergence on the left is cancelled, algebraically, by the 
weight-length-shuffle determination ( |184[ ) of the first term on the right. 

Solving the 36 equations for the 36 complex unknowns, one finds that precisely 5 are 
left undetermined by the pair of shuffle algebras. All 5 remaining real parts are then 
determined by the complex conjugation ( |130| ). Of the 5 remaining imaginary parts, one 
is determined by the generalized duality (|132|) and a second by the generalized duplica- 
tion ( |134| ). Of the 3 that now remain, 2 cannot possibly be determined by the pair of 



shuffle algebras, since ( |166|J169|) entail primitives, for which one requires the analysis of 
section 8.4.3 of 0. 

Finally, just one of 126 reductions achieved empirically awaits rigorous proof. It may 
be taken as ( |170| ). One should not be surprised that the algebra above has, as yet, failed to 
bring this last sheep into the fold of the proven. We know from the example of ( |129| ) that 
there are word transformations which increase the depth of sums. Hence it is eminently 
possible that a relation between double sums requires for its proof the application of the 
machinery of sections 5.1-5.3 to the triple sums to which it may be promoted. 

Here we comment on the reliability of the algebraically unproven and numerically very 
secure result in fll70| ), discovered by applying pslq to 30-digit data. Let us suppose that 
PSLQ is merely a "black box" and ignore its diagnostic messages (which, incidentally, say 
that 30 digits are not needed to be confident of the result). This black box has made a 



prediction, namely that if we compute the 31st digit of the left-hand side of ( |170| ) it will 
agree with that from the right. That would be a weak test, since the probability of a 
bogus reduction agreeing with the new calculation would be 1/10. But now we compute 
the left-hand side to 100-digit precision and find that it agrees with the 70 new digits 
predicted by the right-hand side. If one applied this methodology to reduce 10 70 such 
sums, one would be unlucky to encounter more than one mistake attributable to PSLQ. 



6.2.4 Generalized parity conjecture 

We are led, by the previous results, to suggest a generalization of the parity conjecture for 
EZ sums. In the EZ case, the parity conjecture supposes^ that no new primitive appears 
at depth k and weight n > 1 if the number of f2's, i.e. n — k, is odd. So far our experience 
is consistent with the possibility that the real part of a sum in A\ reduces when n — k 
is odd, and the imaginary part when n — k is even. We refer to this as the generalized 
parity conjecture. The four sectors investigated thus far motivate it, as follows. 

1. At n = 2 and k — 1, the real parts in ( |142|) reduce to tt 2 , while the imaginary parts 
give the primitive Cl2(7r/3). 

8 This is proven for weights n < 8, by massive use of computer algebra. 
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2. At n = 3 and k = 1, the imaginary parts in ( |155f) reduce to 7r 3 , while the real parts 
give the primitive C(3). 

3. At n = 2 and k = 2, the imaginary parts in Table 3 involve no new primitive, while 
the real parts contain Li 2 (j). 

4. At n = 3 and k = 2, the algebra of section 6.2.3 reduces all the real parts, while the 
imaginary parts ( |166j , |169| ) contain the primitives 3Li3(z/y / 3) and ^sLi 3 (X/2). 



Thus at n = k = 3 we expect no new primitive in the imaginary parts and eagerly await 
the appearance of 3^ Li 3 (i/ a/3) and KLi 3 (A/2) as new primitives in the real parts. 



6.2.5 Reductions to exponentially convergent single sums 

It is particularly convenient that the primitive words of A\ at weight n = 3 and depth 
k = 2 are determined by the exponentially convergent single sums ( |158 , 159 ). This finding 



relates strongly to the work on polylogarithmic ladders in ||], where we found transfor- 
mations of EZ sums to polylogarithms inside the unit circle, whose arguments satisfy 
z s = 2~ 4fc with k G {1,3,5}. That work yielding the ten millionth hexadecimal digits 
of C(3) and C(5), just as the simpler results for {it, tt 2 , log2, log 2 2} in |12[ led to rapid 



computation at weights n < 2. We proved that Catalan's constant is also in the class 
SC*(2) of numbers with a <ith binary digit computable in time = 0(d\og T d) and space 
= 0(log 5 d). For the base-2 cases of @ one has @ T < 3 and S = 1. 

Combining the previous and present work, we see a more general pattern of rapid 
computability emerging, with transformations of words in A\ to polylogarithms with 
arguments satisfying z 24 = 2~ J 3 _fe . We were alerted to this possibility at weight n = 2, 
where transformations of dilogarithms [5^] yield 



U 2 (i/V3) = -|(21og2-log3) 2 -iLi 2 (i)-^7rlog3 + f Cl 2 (vr/3) (186) 
Li 2 (A/2) = ^vr 2 + lLi 2 (i)-ivrlog2 + f C1 2 (tt/3) (187) 

whose right-hand sides comprise two- letter words in ^4^- In [0 we showed how to express 
7rlog2 in terms of single sums with z 8 = 2~ 4k . Combining that work with the imaginary 
part of (|187|), we derive the remarkable relation 



= 3 {5 Li 2 (A) - 6 Li 2 (A/2) + 12 Li 2 (i/2) - 16 Li 2 (l±*) + 8 Li 2 (±±i) } (188) 

with only Cl 2 (7r/3) := SLi 2 (A) appearing on the unit circle, while the other terms have 
z 24 = 2~ l2k with k G {1,2,3}. 



Even simpler than (|188|) is the relation of Cl 2 (7r/3) to exponentially convergent dilog- 
arithms with z 12 = 3~ 6 , which is 

= 3 |Li 2 (A) + 6 Li 2 (j\/3 \ - 6 Li 2 (i/V$) | (189) 

yielding the SC*(3) result ( p!7[ ) for the two- loop constant ([T5|). 
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These findings suggest a strategy for identifying primitives in A\. to seek integer 
relations between n-letter words and IA n (z) with z 24 = 2~ J 3~ fc , and more especially z 24 = 
2~i2fe or z i2 _ 3~6fc_ g ere? f or example, are three significant findings at weight n = 3: 

Li 3 (|) = |7r 2 log2-f|C(3) + 8KLi 3 (A/2) (190) 

KLi^^A/?) = - T | i7 r 2 log3 + ^log 3 3 + i|C(3) (191) 

^Li 3 ^/\/3j = -J^ 3 + ^V3+f3Li 3 (i/V3) (192) 

They show that Li 3 (|) and Li 3 (^/A/3) give words of A\, since each reduces to ( |160| ). 

We surmise that the link between DEZ sums and polylogs inside the unit circle is 
similar to that for EZ sums in || : at low weights we find relations between polylogs with 
z 24 = 2~ J 3~ A: ; at higher weights we expect the relations to peter out, with more DEZ 
primitives coming from this source. For the present, there is a strong indication that we 
already have enough primitives on board to complete the investigation at weight n = 3. 



6.2.6 Three-letter triple sums 



Since we are primarily concerned with identifying primitives, we may ignore the product 
terms in the shuffles algebras, and also the final term in ( |118| ), which has lesser depth. 



Thus the algebraic problem for three-letter triple sums is to solve 

) ~ (193) 

((uJpLOqUr) + C(uJqU p UJ r ) + ((uJqUJrUJp) ~ (194) 

where ~ means modulo products and terms of lesser depth. The weight-length shuf- 
fle ( |194j ) is solved by choosing a Lyndon basis. This reduces the 216 words to 70. Solving 
the 216 depth- length shuffles ( |193| ) for the 70 complex unknowns, we find that precisely 
2 are undetermined. Then complex conjugation determines the imaginary parts of the 
remaining two, thereby confirming the prediction of the generalized duality conjecture. 
No further algebraic identity in section 5 determines the two real parts, which are hence 
candidate primitives at n = k = 3. 

To see if the two undetermined real parts are indeed primitive, we may study the pair 
of integrals with a 2 = 1 in 

X a :=\^dy\og{l + ay + y 2 )±\og 2 (^) (195) 

which are easily evaluated to 100 digits, as they have good behaviour at the endpoints. 
Then pslq finds that 

3?C(^ 2 ) = X 1 = ^7r 2 log2-|log 3 2-|C(3)+KLi 3 (A/2) (196) 
KCKV) = X_ 1 = T | i vr 2 (41og2-log3)-^(121og 2 2-log 2 3)log3 

-^C(3)-X!-^Li 3 (i/V3) (197) 



confirming that ( |160| ) is a viable basis at n = 3. By rigorous computer algebra we 
have shown that there are no more than 5 new constants entailed by double and triple 
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sums of weight n = 3. The 5 empirical reductions fll66| , 169 , 170| , 196 , 197 ) then relate 
the algebraically undetermined terms to £(3) and the real and imaginary parts of the 
two exponentially convergent complex trilogarithms in (|160|) . Moreover, pslq found no 
integer relation within this set of constants, which is the best that one can hope to say 
about primitives, in the foreseeable future. 

Thus we are done with three-letter words. There is one real primitive, ((3) at depth 
k — 1, and two complex primitives, Li 3 (z/\/3) and Li 2 (A/2), whose imaginary parts enter 
at k — 2, with the real parts postponed until k = 3, in line with the generalized duality 
conjecture. Neither the Deligne sub-alphabet, A\ := {Q, u> , u^, nor the Euler-Zagier sub- 
alphabet, A3 := {Q,uj ,uj 3 }, entails the complex pair. The present analysis, combined 
with || , results in the remarkable finding that all 800 real and imaginary parts of words in 
A\ with less than 4 letters have finite parts reducible to SC* constants and their products. 



6.3 Four-letter words in field theory 

In section 6.2, life was comparatively easy: all words with 3 letters were expressible as 
single integrals of the form Jq 1 dy A(y) B(y) dC(y)/dy, where A, B, C were simple logs, or 
their imaginary parts: arctangents. Now, with 4-letter words, we encounter less pleasant 
things, since one of the trio A,B,C becomes a dilog. If it is C, things are not too bad. 
If it is A or B, as occurs more often, the time required to compute the word increases. 
Moreover we need the results to higher precision than before, because the size of the basis 
has grown, due to the proliferation of product terms. Three things help us to surmount 
these computational problems. 

First, the pair of shuffles algebras greatly reduces the number of words that we need to 
study, allowing us to choose the undetermined terms as those most easily calculated. Sec- 
ondly, maple has a convenient implementation of high-precision one-dimensional quadra- 
ture, enabling fast evaluation when the endpoint behaviours are benign. Finally, PSLQ 
revealed yet more features of the extraordinary relation between primitives of the shuffle 
algebras, Feynman diagrams, and super-fast computability in bases 2 and 3. 



6.3.1 Rapid computability of the Clausen-Deligne primitive 

For weight-4 sums with depth k — 1, we obtain 4-letter words with 3 fi's: 

for p = 



C(n 3 cu p ) = Li 4 (A p ) = < 



' u vr 4 + zC1 4 (tt/3) for p = 1 
, )(30 7r 4 + |Cl 4 (7r/3) forp = 2 

720 



19440 
13 ^4 



(198) 



for p 



where Cl4(7r/3) := J2 n >o sin(7rn/3)/n 4 is presumed primitive. Its occurrences at p = 1 
and p = 2 are related by 



Cl 2m (7r/3) 



)l-2m 



Cl 2m (27r/3) 



(199) 



which was used in ((L42|) and in the definition (|15|) of the constant S2 in the finite part of 
the 2- loop vacuum diagram with 3 masses [^, |TI], [□]] . 
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In general, the real parts at depth k = 1 and even weight n = 2m come from [22] 

2{-A) m - l B 2m forp = 

( 2 l-2m _ 1)(3l-2m _ ^ {-A)" 1 ' 1 B 2m for p = 1 



7T-'" (S 1 " 2 ™ - l)(-4) m - 1 J B 2m forp = 2 ^ 20 °' ) 

(2- 2 ™-2)(-4) m - 1 J B 2m forp = 3 

with the Bernoulli number B4 = — ^ giving the real parts in ( 198 ). In the case of 
C,{Vt 2m ujp) one encounters ((2m + 1) in the real parts and 7r 2m+1 in the imaginary parts. 
Hence, at each weight n > 1, a single depth-1 constant is presumed to be primitive: at 
odd n it comes from Riemann; at even n from Clausen. Euler-Zagier sums are real and so 
do not involve the Clausen sum, which is however entailed by the Deligne sub-alphabet, 
at depth k = 1 and, by generalized duality, at maximum depth. Hence we characterize 

- 9f CH^o 2 - 1 ) = 9 C(fi 2m ~ V) = Cl 2m (7r/3) := £ (201) 

as a Clausen-Deligne primitive of depth k — 1 and weight n = 2m. 

The next objective is to find a computationally more convenient representation of 
04(71-/3), as the imaginary part of an exponentially convergent quadrilogarithm, with 
z 24 = 2~ J 3~ fc , inside the unit circle. We then expect to encounter the real parts of such 
quadrilogarithms as primitives in the Feynman sector, with weight n = 4 and depth k = 2. 
At weight n = 2, the dilogs ( |186| , |187| ) serve a similar purpose, since each generates both 
two-letter primitives. We do not expect the n = 4 case to show so much redundancy. 

As an intermediate step, we can use the representation 

with finite sums /i r (n) := I] 2 ™^ 1 s~ r , obtained by expanding in y(l — y) and taking the 
third differential of an Euler Beta function. Updating all values in a single exponentially 
convergent loop, we obtained a 1000-digit value in a few seconds. 

As a target we chose the 5 independent imaginary parts of Li 4 (z) with z 12 = 3 -6 , since 
2 such arguments figured in (|189|) at weight n = 2, and in (|192|) at n = 3. Working at 
150-digit precision, pslq found that 

v / l\ n ( 9 15 18 5 1 1 

n V \~27j { (6n + l) 4 ~ (6n + 2) 4 ~ (6n + 3) 4 ~ (6n + 4) 4 + (6n + 5) 4 J 

= V? (f Cl4 ( 7r / 3 ) " f^ 3 lo S 3 + S w lo S 3 3 ) ( 203 ) 
with an element of SC*(3), on the left, expressed in terms of the imaginary parts of 4-letter 
words, up to the predictable factor of sin(7r/3), already observed in the imaginary parts 
of ( |158|J159|) at weight n — 3, and in definition (|15[) at n — 2. We have checked that (|203|) 
correctly predicts the next 850 digits of 04(71-/3), obtained from the proven result ( |202| ). 

It is instructive to compare ( [203]) with the corresponding result at weight n = 2: 

v / l\ n ( 9 15 18 5 1 ^ 



n>0 



27/ [(6n+l) 2 (6n + 2) 2 (6n + 3) 2 (6n + 4) 2 (6n + 5) 2 J 

3 7rlog3 (204) 



33 



which shows that sin(7r/3) 7rlog3 e SC*(3). Already we see a realization of the hope 
that the relation between SC*(3) and A\ will help us to identify primitives: the integer 
coefficients on the left of (|204|) that give merely a product term at n = 2 yield a primitive 
in ( |203| ) at n = 4. This observation proves extremely useful in section 6.3.3 

Having thus found the connection of quadrilogs with z 12 = 3 -6 to the 4-letter Clausen- 
Deligne primitive at depth k — 1, we move on to study the real parts at k — 2, where 
Feynman diagrams have already generated two presumed primitives: f/ 3>1 and V^i in the 
basis (|116|) for the vacuum diagrams of Table 2. 



6.3.2 Feynman rules OK: no more depth-2 primitives 

There are 108 four-letter depth-2 words. The 36 of type ((u p Q 2 uj r ) may be eliminated 
by the weight-length shuffles of ((u> p ) with ((Q 2 u r ); the 36 of type ((flu p flu r ) by the 
corresponding depth-length shuffles. The remaining 36 of type ((Q 2 u p uj r ) are reduced to 
7 by the weight-length and depth- length shuffles of ((Qlv p ) with ((Qu r ). The 7 remaining 
imaginary parts are then determined by the complex conjugation ( |130| ), leaving only real 
parts, consistent with the generalized parity conjecture. Of these 7 real parts, one is 
determined by the generalized duplication ( |134j ) and another by the triplication ( p.35| ). Of 
the 5 real parts that now remain, we account 

( 1 \m+n 

U 3A := £ ^-L = ((flWo) (205) 

m>n>0 m 71 

_ (-l) m cos(2vm/3) , . 

^3,1 := E — r — = ^C(^3UJi) (206) 

m^>0 771 71 

as Feynman primitives, since we encountered them in the calculation of vacuum diagrams, 
computed them to 1000 digits, and found no reduction to products, using pslq. 

There thus remain 3 real parts to consider. They are undetermined by the algebraic 
methods used so far, but need not be new primitives. In section 6.2.3 we noted that there 
are relations beyond those derivable from shuffles, complex conjugation, duplication and 
triplication, attributable to the interplay of these with word transformations ( |127|jl32| ). 
since the latter change the number of Q's and hence the depth. It follows that at a given 
depth there can be reductions that are apparent only after one applies the full algebraic 
machinery to all depths. This is where PSLQ is of such enormous benefit: rather than 
having to solve, simultaneously, for 2 x 7 4 = 4802 unknowns, in the real and imaginary 
and imaginary parts of all 4-letter words, we break up the problem by depth, relying on 
the lattice algorithm to find the small percentage (in this case ^ < 1.5%) of reductions 
that we miss by staying at fixed depth. Since we miss so little by staying in a fixed sector, 
the small amount of numerical investigation can focus on the most easily computable of 
those sums which partial (though already very powerful) algebra has left undetermined. 

Accordingly, we evaluated to 100-digit precision the 4 conveniently well-behaved, dilog- 
free, quadrilogarithmic integrals with o 2 = t 2 = 1 in 

Y a>T :=\ £ dy log 2 (y)\og(l- ay) -^\og(l + ry + y 2 ) (207) 

with the duplication relation Y 1;1 = 4 J2ar Y a , T providing a good check of accuracy. The 
algebraic reductions of 211 of the 216 real and imaginary parts then prove that the question 
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whether the 2 Feynman primitives suffice, at depth k = 2, is identical to asking whether 
3 instance of (|207| ) are reducible to ( |205| , [206| ) and/or product terms. According to pslq, 
Feynman rules OK, with 

&C(tfW*>) = r 1)1 = ^7r 4 -§ 01^/3) + f^,i (208) 
^C(^V^o) = lV 1 = 3^7r 4 (209) 
KC(^W 3 ) = F-i,i = -2i27r 4 + |C(3)log2-|f/3,i (210) 

»C(nV 3 ) = = -111,1-11,-1-^-1,1 (2ii) 

found easily at 50-digit precision, and then tested stringently by the next 50 digits. 

Note that the simple reduction (|209|) is entirely as expected: the word is in the Deligne 
sub-alphabet, where the generalized duality ( |132| ) immediately promotes it to depth k = 3. 
We have every confidence that its reducibility is provable by application of algebra at 



depths k > 2, just as (|165|) , obtained empirically at depth k = 2 in section 6.2.1, was 
later proven by algebra at depth k = 3 in section 6.2.2. However, we are content to 
have completed the analysis of over 1 000 real and imaginary parts relevant to the 3-loop 
Feynman sector, by restricting attention to weights n < 4, and to depths k < 2 at weight 
n = 4. We beg the reader's pardon for postponing analysis of over 4 000 cases, with n = 4 
and k > 3, to a later occasion, when quantum field theory is not the pressing issue. 

In conclusion: we are done with enumeration of primitive 4-letter words at depth 
k = 2. There are precisely two. Feynman found them both, way back in Table 2. Neither 
is in the Deligne sub-alphabet. Neither occurs with 3 masses. One is in the Euler-Zagier 
sub-alphabet and occurs only when the number of masses is even; the other is not in the 
Euler-Zagier sub-alphabet and occurs only with 5 masses. 

6.3.3 Number and topology in the Feynman sector 

In the Feynman sector, with weight n = 4 and depth k = 2, two questions remain: 

Ql: Are the quadrilogarithms from all 10 Feynman diagrams reducible to elements of 
SC*(3) U SC*(2) and their products? 

Q2: If so, how does the base of super-fast computation, b = 2 or b = 3, relate to the way 
in which mass colours the Feynman tetrahedron? 



The first issue to address is the status of 7r 4 . We know from that ir A is the square 
of an SC*(2) constant, and very recently learnt in || that 7r 4 e SC*(2), which is a highly 
non-trivial result, since there is no indication that multiplication is allowed in SC*. Indeed 
the polylogarithmic ladders of || suggest that 7r 6 is not in SC*(2). Now we add 



72 

(1271 + 6) 1 
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to the jigsaw puzzle. It shows that it 4 is the square of a constant in SC*(3). So is 
Cl^(7r/3), since (0) shows that v/3C1 2 (tt/3) G SC*(3). From we conclude that 

G SC*(2). Hence we need only investigate the 5- mass primitive Vjj^. 

There is a clear SC*(3) candidate for a relation to V31. Noting that an increase in 
weight from n = 2 in ( |204| ) to n = 4 in ( |203| ) promoted a product to a primitive, we 
increased the weight in ( 212 ) and obtained V^x in association with products, with 

ivr 4 + ffvr 2 log 2 3 + !f (01^/3) - 3 V 3>1 



48 1 16 to 1 2 

y, / 1 \ n | 243 405 81 27 

n>0 ^729/ [ 



(12n+l) 4 (12n + 2) 4 (12n + 4) 4 (12n + 5) 4 
72 9 9 5 1 



(12n + 6) 4 (12n + 7) 4 (12n + 8) 4 (12n+10) 4 (12n + ll) z 



(213) 



checked to 1000 digits, using the value for V^x from section 6.2.1. Since 7r 4 and Cl2(vr/3) 
are squares of SC*(3) constants, only the status of ir 2 log 2 3 in ( |213[ ) remains an issue. It 
is clear that 

/ 1 f 729 81 81 

^§3 = E ™ 77^T + 77— T^7 + 



n>0 



V729/ [(6n+l) (6n + 2) (6n + 3) 

9 9 1 1 

+ (6ra + 4) + (6n + 5) + (6n + 6) J 

is an SC*(3) constant, like the weight-1 constant a/37t of massive Feynman diagrams f37j| . 
We now wind ( pi4| ) up to weight n = 2, where the corresponding relation is 

1 , 20 , /nN / 1 f 729 81 81 

±log 2 3 = C(2)-£ — 77—^ + 77—^ + 



n>0 



V729/ [(6n + l) 2 (6n + 2) 2 (6n + 3) 2 

9 9 1 1 

+ (6n + 4) 2 + (6n + 5) 2 + (6n + 6) 2 J ( } 

which shows that log 2 3 G SC*(3), since (0^) gives C(2) = |vr 2 G SC*(3). Thus (|2T3| ) 
shows that the second Feynman primitive, V^i, consists of an SC*(3) constant and prod- 
ucts of pairs of SC*(3) constants. So therefore does the quadrilogarithm in the 5-mass 
diagram. Thus the answers to the questions above are as follows. 



Al: Remarkably, all 10 quadrilogarithms reduce to elements of SC*(3) U SC*(2) and 
their products. It follows that evaluation of their decimal values is fairly (though 
not super) fast. We obtained 10 000 decimal places, at a cost of 25 seconds/diagram, 
on a 333 MHz workstation. 

A2: The quadrilogarithms in {Vi,V 2 a,V3t} are multiples of 7r 4 G SC*(2), which is also 
the square of n 2 G SC*(3). Those in {V^V^v} are constants in SC*(2); those in 
{Vss, V3L, V^a} are sums of squares of constants in SC*(3). The 5-mass case reduces 
to SC*(3) constants and their products; the 6- mass case to an SC*(2) constant and 
the square of an SC*(3) constant. The primitive U^x G SC*(2) does not occur when 
the number of masses is odd; the other primitive, V^i, is unique to the 5-mass case. 
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7 Conclusions 



We have come a long way, from the innocent idea of identifying analytically the unknown 
quadrilogarithm in the 3-loop QCD corrections to the rho-parameter 0, ||] of the standard 
model, to an investigation of the shuffle algebras of polylogarithms of the sixth root of 
unity, revealing a very small number of primitives that take values from constants in 
SC*(3)USC*(2) and their products. Along that route, it was quantum field theory that led 
us, by wonderfully simple results of decidedly non-trivial calculations, still only partially 
elucidated by painstaking analysis. The following features seem the most remarkable. 

First, who would have imagined that three- loop Feynman diagrams with 3 masses 
are simpler than those with 2 masses? With 2 masses, Table 2 shows that we already 
engage with a primitive four-letter word of the sub-alphabet of Euler-Zagier sums; yet 
with 3 masses we merely pick up the squares of constants already encountered as two- 
letter words, at two loops. The origin of the simplicity of 3-mass diagrams appears to 
be the superiority of the primitive sixth root of unity over the primitive square root: 
the Deligne sub-alphabet has a simpler generalization ( |132| ) of the Zagier duality ( |131| ) 
than that in (|127|) for the Euler-Zagier sub-alphabet. It seems that the star-triangle-line 
relation flgU| ) is an expression of this: real parts of 4-letter Deligne words entail only the 
products 7r 4 and Cl^vr/S) at depth k = 2; these are the only constants that appear with 
3 masses, which is why a relation exists. The result, 3 stars = triangle + 2 lines, cries out 
for an explanation of why 3-mass diagrams entail only the Deligne sub-alphabet and why 
the integers in the ensuing relation have such simple values. Such simplicity is far from 



apparent in (|68|), proven by a dispersion relation. 

Next, who would have imagined that with 4 masses no new constant is entailed? 
Table 2 shows that it is the topology of the diagram that determines the previously known 
case to which one reverts. When the two massless lines are adjacent one reverts to the 
simplicity of 3-mass cases, entailing reducible words from the Deligne sub-alphabet; when 
they are non-adjacent one reverts to the primitive four-letter word of the Euler-Zagier 
sub- alphabet, as found in the 2- mass case when the massive line were non-adjacent, with 
now a double measure, for double the number of masses. 

Then, who would have imagined that while 5 = 3 + 2 masses entails a new primitive, 
6 = 3x2 masses is far simpler? The totally massive case yields only the Euler-Zagier 
constants £(4) and U^i, already found with 2 masses, and the square of the Clausen- 
Deligne constant Cl2(vr/3), already seen with 3 masses, at merely two loops. This last 
finding shows that the current state of the art lags far behind the simplicity of the results. 
Until this work, it had proven too arduous an undertaking to identify even the simple 

sum of squares (Cl2(7r/3)) 2 + (|C(2)) = C(3) — \D% in the 3-mass quadrilogarithm that 
provides the value for D 3 := V^l in the 3-loop QCD corrections to the rho-parameter ||, 
Despite the much more formidable obstacle posed by the Cutkosky rules for the 6-mass 
case, double integration by nag and an integer-relation search by pslq finally revealed the 
totally massive beast to be a beauty. It was scarcely to be expected that the convolution 
of a dilogarithm with an elliptic integral of the third kind would yield merely the unit 
coefficients of the integer relation ( |115| ), reducing the 6-mass case to those with 4 or less. 



Now one sees why Leo Avdeev's thorough analysis of integration by parts || found no 
route from 6 masses to 5, nor from 5 to 4. There can be none, since the 5-mass case is 
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the only one to entail the second Feynman primitive. 

Finally, who would have imagined that the primitives left undetermined by the shuffle 
algebras would entail only the classes SC*(2) and SC*(3) of super-fast computability? 
Diagrams with 2 masses choose SC*(2); those with 3 choose SC*(3); that with 6 = 3x2 
masses chooses SC*(3) U SC*(2). One clearly needs to decode this mapping from the 
10 colourings of a tetrahedron, by mass, to the base of super-fast computation of the 
constants that field theory assigns to those colourings. 

Truly the adage "Out of field theory always something new" has been confirmed. 
I hope that I have spelt out enough of the very recently discovered detail to convince 
mathematical colleagues that we physicists sit on top of a structural gold mine, and fellow 
physicists that we have much to learn about the relation between the free investigations of 
mathematics and the tightly constrained structure of perturbative quantum field theory 
in four-dimensional spacetime. 
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